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ABSTRACT 



' Context. Generation and dissipation of magnetic fields is a fundamental physical process on the Sun. In comparison to flux emergence 

I and the initial stages of sunspot formation, the demise of sunspots still lacks a comprehensive description. 

• , Aims. The evolution of sunspots is most commonly discussed in terms of their intensity and magnetic field. Here, we present additional 

• information regarding the three-dimensional flow field in the vicinity of sunspots towards the end of their existence. 

Methods. We present a subset of multi-wavelengths observations obtained with the Japanese Hinode mission, the Solar Dynamics 
JL I Observatory (SDO), and the Vacuum Tower Telescope (VTT) at Obsen'atorio del Teide, Tenerife, Spain during the time period from 

^ , 2010 November 18-23. Horizontal proper motions were derived from G-band and CaiiH images, whereas line-of-sight velocities 

were extracted from VTT Echelle Ho- /1656.28 nm spectra and Fei /^630.25 nm spectral data of the Hinode/Spectro-Polarimeter, 
which also provided three-dimensional magnetic field information. The Helioseismic and Magnetic Imager on board SDO provided 
continuum images and line-of-sight magnetograms as context for the high-resolution observations for the entire disk passage of the 
active region. 

Results. We have performed a quantitative study of photospheric and chromospheric flow fields in and around decaying sunspots. 
^ ' In one of the trailing sunspots of active region NOAA 11126, we observed moat flow and moving magnetic features (MMFs), even 

^\ ' after its penumbra had decayed. We also noticed a superpenumbral structure around this pore. MMFs follow well-defined, radial 

00 , paths from the spot all the way to the border of a supergranular cell surrounding the spot. In contrast, flux emergence near the other 

l/^ . sunspot prevented it from establishing such well ordered flow patterns, which could even be observed around a tiny pore of just 

2 Mm diameter After the disappearance of the sunspots/pores a coherent patch of abnormal granulation remained at their location, 
which was characterized by more uniform horizontal proper motions, low divergence values, and diminished photospheric Doppler 
velocities. This region, thus, differs significantly from granulation and other areas covered by G-band bright points. We conclude that 
this peculiar flow pattern is a signature of sunspot decay and the dispersal of magnetic flux. 

Key words. Sun: chromosphere - Sun: photosphere - Sun: surface magnetism - Sun: sunspots - Techniques: image processing - 
Methods: data analysis 



H ' 1. Introduction entire spot, was hindered by continuous flux emergence between 

. . 1 the spots of the bipolar region. 

Sunspots are a thought-provoking aspect of solar activity be- Quite the opposite, the decay of a sunspot is a slow process, 
cause of the close interaction between plasma motions and mag- Decay rates for stable leading sunspots and irregular follower 
netic fields. Recent progress in MHD simulations (e.g., Rempel, spots are different (Martinez Pillet, 2002). A number of decay 
2011) provide a comprehensive framework for the interpre- laws were proposed such as a linear decay law described by 
tation of high-resolution sunspot observations. The formation Bumba (1963) and a parabolic decay law proposed by Petrovay 
of a penumbra around a sunspot is a rapid phenomenon, i.e., & van Driel-Gesztelyi (1997). Martinez Pillet (2002) critically 
within a few hours a sunspot can develop a penumbra (Leka reviews various diffusion models, concludes that they explain 
& Skumanich, 1998; Yang et al., 2003), which is intimately well how flux is spread over larger areas while the spot is de- 
linked to more inclined magnetic field lines and the onset of the caying, but they fail to satisfactorily describe the flux removal 
Evershed flow. Schlichenmaier et al. (2010) observed the growth process. The initial stages of sunspot decay, i.e., when the spot 
of a penumbra where only the newly formed penumbra con- looses its penumbra, are exemplary described in Bellot Rubio 
tributed to the increase in spot size while the umbra remained et al. (2008), who discovered finger-like structures, which are 
stable. The formation of a penumbra, which would surround the neither related to penumbral filaments nor the Evershed flow. 
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Fig. 1. Calibrated Hinode G-band images showing the decay of active region NOAA 111 26 observed during the period from 2010 
November 18-23 (from top-left to bottom-right). The FOV is 111" x 111". The annotation of the axes refers to heliocentric co- 
ordinates given in seconds of arc. Brighter areas delineated by vertical black lines correspond to regions, which were covered by 
spectral scans obtained with the VTT Echelle spectrograph. The intensity scale to the right applies to these regions, while areas 
not covered by Echelle data are displayed with an offset of O.I/q. On 2010 November 23, only few Hinode G-band images were 
available. The limb darkening has been subtracted from the G-band images, which were then normalized so that the mean of the 
quiet Sun intensity distribution corresponds to unity. 



These features might be penumbral field lines rising to the chro- 
mosphere, thus contributing to the decay of the sunspot penum- 
bra. When a sunspot looses its penumbra, its decay reaches a 
critical point. Magnetic field lines become more vertical and 
convective motions in its vicinity begin to change. These ideas 
of a critical inclination angle and convective motions were put 
together by Rucklidge et al. (1995), who explains in a simple 
model why small sunspots can have a penumbra while larger 
pores do not possess one. 

The moat flow is a large-scale flow pattern commonly ob- 
served around sunspots (Meyer et al., 1974). However, flux re- 
moval and dispersal can only be understood in the context of the 
moat flow's fine structure. Moving magnetic features (MMFs) 
play a major role in the flux dispersal process and they are only 
associated with decaying sunspots (Harvey & Harvey, 1973). 
The total flux carried by MMFs is several times larger than the 
flux contained within the sunspot itself. Thus, the polarity of 
MMFs has to be considered for a balance of the net flux. MMFs 
move radially outward with a velocity of 1 km s"' before they 
reach and dissolve within the network, i.e., at the boundaries of 
the supergranular cell containing the sunspot. Zuccarello et al. 
(2009) showed evidence that MMFs and moat flow are present 
even in the vicinity of pores, i.e., in the absence of penumbral 
filaments and Evershed flow (cf., Cabrera Solana et al., 2006). 
Deng et al. (2007) also detected a persistent moat flow after the 
penumbra around spot disappeared leaving only a pore. Even 



though the moat flow might not be closely tied to the Evershed 
flow, the sub-photospheric interaction of magnetic field lines and 
flows can still produce the observed flow patterns. 

Verma & Denker (201 1) described a local correlation track- 
ing (LCT) method to measure horizontal flows based on Hinode 
G-band images. In this study, we perform a case study, where 
we put such horizontal flow fields in the context of other pho- 
tospheric and chromospheric data. In particular, we are inter- 
ested in the final stages of sunspot decay. In Sect. 2, we present a 
subset of multi-wavelengths observations, which were obtained 
within the scope of Hinode Operation Plan (HOP) 0176. The 
temporal evolution of active region NOAA 11126 in terms of 
intensity, morphology, and flow as well as magnetic fields is de- 
scribed in Sect. 3 and discussed in Sect. 4. 



2. Observations 

The disk passage of active region NOAA 1 1 126 started on 2010 
November 12 and ended on November 24. NOAA 11126 was 
classified as a y6-region, while it crossed the solar disk. No major 
flaring was associated with the region. Only a few B-class events 
were reported on 2010 November 15. As part of HOP 0176 
''High- resolution multi-wavelength study of small-scale jets on 
the solar disk", we observed the decay of two small follower 
sunspots in the active region for five days from 2010 November 
18-22. A time-series of Hinode G-band images is shown in 
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Table 1. Observing characteristics and physical parameters 
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The parameters in the first column of the table refer to 
heliographic latitude B. heliographic longitude L, cosine 
of the heliocentric angle start of the G-band (GB), 
spectro-polarimeter (SP), and Echelle spectrograph (ES) 
observing sequences fo, spot area A, horizontal velocity v, 
mean divergence |Vii, number of magnetic elements Wmas- 
The indices refer to G-band bright points (bp), granula- 
tion (gran), magnetic elements (mag), the two sunspots 
(Spot A and Spot B), and the four-megameter wide an- 
nuli around both spots (Ring A and Ring B). If two rows 
are given for a physical parameter then the top row refers 
to G-band data, whereas the bottom row was derived from 
Ca n H data. If present, the standard deviation refers to the 
variation of the physical parameters within the specified 
regions rather than to any formal en'or 



Fig. 1, where we labeled the northern and southern spots with 
A and B, respectively. In the following, we will simply refer to 
these magnetic features as spots, even if the proper classifica- 
tion should be pores, i.e., sunspots lacking a penumbra. Spectral 
scans with different field-of-views (FOVs) and cadences were 
observed for three hours every day. Since we are only focusing 
on the general properties of sunspot decay, we chose the first 
scan on the given day, which covered the largest FOV. The set- 
tings for the Vacuum Tower Telescope (VTT) Echelle data and 
data of Hindoe/Spectro-Polarimeter- were chosen as to obtain 
the best spatial and spectral match. The general observing char- 
acteristics are listed in Tab. 1 . 



2.1. SDO/HMI full-disk images 

The discussion of the temporal evolution and morphology of ac- 
tive region NOAA 1 11 26 is based on full-disk images and line- 
of-sight (LOS) magnetograms obtained with the Helioseismic 
and Magnetic Imager (HMI, Schou et al., 2010; Couvidat et al., 
2011; Wachter et al., 2011) on board the Solar Dynamics 
Observatory (SDO). Since the Hinode FOV is too small to cover 
the entire active region and to give an overview of the magnetic 
field topology, we show in Fig. 2 the limb-darkening corrected 
HMI continuum image and magnetogram for November 18. 

We selected from the SDO/HMI database one im- 
age/magnetogram with 4096 x 4096 pixels every 15 min for the 
period from 2010 November 13-23, i.e., a total of 1056 full- 
disk images. The image scale is about 0.5" pixel so that even 
finer details of penumbra and umbra can be captured. The aver- 
age limb-darkening function was computed for this time interval 
and subtracted from the full-disk images to yield contrast en- 
hanced images (see e.g., Denker et al., 1999), which can then be 



used for feature identification. The photometric temporal evolu- 
tion for the entire active region as shown in Fig. 2 is depicted 
in Fig. 3. The corresponding changes of the magnetic flux are 
shown in Fig. 4. Since HMI data cover the whole solar disk, it is 
straightforward to compute the heliocentric angle yU on a pixel- 
by-pixel basis. Thus, a geometrical correction is applied to the 
measured areas and flux values of Figs. 3 and 4, which are dis- 
cussed in detail in Sects. 3.1 and 3.5, respectively. 

2.2. Hinode G-band and Ca ii H images 

We applied LCT (for details see Verma & Denker, 2011) to 
image sequences captured in G-band /1430.5 nm and Can H 
/1396.8 nm to compare horizontal flows in the photosphere and 
chromosphere. Note that the Ca ii H images do not purely repre- 
sent the chromosphere, but contain contributions from both the 
upper photosphere and lower chromosphere. These observations 
were carried out by the Broad-band Filter Imager (BFI) of the 
Solar Optical Telescope (SOT, Tsuneta et al., 2008) on board 
Hinode (Kosugi et al., 2007). Data sequences were captured ev- 
ery day from 09:00 UT to 12:00 UT with an average time ca- 
dence of 120 s (with some jumps in the data sequences). In both 
wavelengths the images are 2 x 2-pixel binned with an image 
scale of 0. 11" pixel"' . Images have a size of 1024 x 1024 pixels 
and a FOV of 111" X 111". 

After basic data calibration, the images were corrected for 
geometrical foreshortening and resampled onto a regular grid of 
80 km X 80 km. The signature of the five-minute oscillation was 
removed from the images by using a three-dimensional Fourier 
filter with a cut-off velocity of 8 km s ' corresponding roughly to 
the photospheric sound speed. For measuring horizontal proper 
motions, we applied the LCT algorithm described in Verma & 
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Fig. 3. Temporal evolution of the area covered by active region 
NOAA 1 1 126 while it crossed the solar disk. The area enclosed 
by umbral cores and pores is displayed with dashed line. Dashed 
dot dot line refers to the penumbra, whereas solid line denotes 
the total area. Some smoothing was applied to the time-series to 
suppress features on temporal scales below one hour. The verti- 
cal gray bars refer to the observing periods of HOP 0176 (2010 
November 18-22). The labels indicate different stages of active 
region evolution, which are explained in Sect. 3.1. 



Fig. 2. Limb-darkening corrected HMl continuum image (top) 
and magnetogram (bottom) of active region NOAA 11126 on 
November 18. The square box in both images shows the FOV 
covered by Hinode/BFl. The axes are labeled in heliographic co- 
ordinates. 



Denker (201 1), which computes cross-correlations over 32 x 32- 
pixel regions with a Gaussian kernel having a FWHM of 15 pix- 
els (1200 km) corresponding to the typical size of a granule. In 
two aspects we deviated from the aforementioned algorithm, the 
time cadence was Af = 120 s and the flow maps were averaged 
over AT = 3 h. 

2.3. Ha Echelle spectra 

The observations in Ho- /1656.28 nm and Fei /1656.92 nm were 
carried out with the VTT Echelle spectrograph. Spectral data 
were acquired with a slit width of 80 fim and an exposure time 
of 300 ms. The image-scale of the spectrograph is 8.99" mm"'. 
We did not use a predisperser, hence, to suppress overlapping 
in spectral orders, we placed a broad-band interference filter di- 
rectly behind the spectrograph slit. The infrared grating with a 
blaze angle of 51.6° and 200 grooves mm"' was used to record 
spectra in the 12* order In this configuration, we achieved a 
dispersion of 0.60 pm pixel"' . The spectra covered a wavelength 
range of 1.2 nm from /1655.9 nm to A651.1 nm. We employed a 
PCO.4000 CCD camera that has a quantum efliciency of about 
30% at Ha. After 2 x 2-pixel binning the images have a size of 
2004 X 1336 pixels. The pixel size of the CCD detector is 9 /im 
X 9 yum. The Echelle data were intended to match the Hinode 
observations, which was mostly achieved except for a few in- 
terruptions because of deteriorating seeing conditions. The two- 
dimensional FOV was scanned with a spatial step of 0.32" and 
200-250 spectra were recorded in a sequence. A sequence of 240 
spatial steps took about 12 min and covered a FOV of 72.0"x 
182.6". 

The common FOVs of the G-band images and Ha Echelle 
spectra are shown in Fig. 1 for each observing day. We first 
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Fig. 4. Temporal evolution of the magnetic flux contained in ac- 
tive region NOAA 11126 while it crossed the solar disk. The 
solid, dashed dot dot, and dashed lines refer to the total, posi- 
tive, and negative flux, respectively. 

matched the image scale of the Echelle data to that of the Hinode 
data. We then aligned Hinode G-band images and continuum im- 
ages derived from Echelle spectra. After this procedure, the he- 
liocentric coordinates for G-band images and Ha spectra differ 
by less than 1" in the periphery of the FOV. Furthermore, the 
Echelle scan direction is not perfectly aligned with the Hinode 
BFI detector. Hence, we computed an offset angle (smaller than 
+2°) for each date and applied it to the spectral data. On 2010 
November 20, the time difference between G-band images and 
Echelle spectra was about five hours. 

2.4. Hinode Spectro-Polarimeter 

The photospheric magnetic topology and evolution of the active 
region were studied using high-resolution spectral data from the 
Hinode/ Spectro-Polarimeter (SP, Ichimoto et al., 2008), which 
uses two magnetically sensitive Fei lines at /1630.15 nm and 
/1630.25 nm and the nearby continuum (Tsuneta et al., 2008) 
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to obtain Stokes IQUV spectral profiles. We used spectral data 
captured in the fast mapping mode with a FOV of 58"x 122" and 
an average scan time of 12 min. On November 18, we also used 
ten continuous high cadence scans with a FOV of 32"x 123" and 
a scan time of 7 min. The dispersion is about 2.155 pm pixel"'. 
The region was scanned with a spatial step of about 0.3" and 
image scale of 0.32" pixel"'. The basic data reductions such as 
subtraction of dark current, flat fielding, polarization, and wave- 
length correction were performed using procedures available in 
SolarSoft (SSW, Bendey & Freeland, 1998; Freeland & Handy, 
1998). 

3. Results 

3. 1 . Photospheric evolution 

Contrast-enhanced HMI full-disk continuum images were used 
to study the evolution and decay of active region NOAA 1 1 126 
during its disk passage. The curves in Fig. 3 correspond to the 
areas of umbral cores/pores, penumbrae, and the sum of both 
types of features. These strong magnetic features are identified 
according to intensity thresholds of 75% and 92%, respectively, 
where the quiet Sun intensity was normalized to unity. Some spa- 
tial smoothing and minimum-size criteria were applied to binary 
masks of identified features using morphological image process- 
ing techniques, thus ensuring unwrinkled boundaries and con- 
tiguous structures. Note that this algorithm only provides a rough 
estimate of the above areas. Some (small) features could be mis- 
classified. In particular, the borders of pores will be classified as 
penumbra. 

The most important stages of the active region evolution are 
labeled in Fig. 3: (a) Initially, two tightly spaced sunspots of pos- 
itive polarity were present early on November 13. The leading 
spot was larger and had a well established penumbra, (b) New 
flux emerged towards the south-east of these spots at 8:30 UT 
on November 14. Numerous (more than ten) magnetic knots and 
pores appeared to the south-west forming a bipolar magnetic re- 
gion, (c) The umbral core of leading sunspot of the new group es- 
tablished at 2:00 UT on November 15. (d) The leading sunspot of 
the new group continuously grew by advecting magnetic knots 
and small pores. The umbral cores/pores occupied the largest 
area at 2:20 UT on November 16. (e) The penumbra of the lead- 
ing spot reached its maximum about nine hours later. At this 
time, the active region NOAA 1 1 126 started its decay phase, (f) 
Some further flux emergence occurred in the trailing part of the 
region at about 16:00 UT on November 17, which strengthened 
spot B and produced thin elongated dark lanes. These typical 
features of flux emergence (see e.g., Strous et al., 1996) were la- 
beled erroneously as penumbrae by the thresholding algorithm. 

HOP 0176 focused on the two trailing spots/pores, of which 
the northern one (spot A) was already decaying, while the south- 
ern one (spot B) had just reached its maximum. This sunspot also 
showed strong indications of rotation, (g) This spot then frag- 
mented into numerous magnetic knots until about 19:30 UT on 
November 19. (h) At this time, the fragments started to converge 
again forming a small sunspot, which reached its maximum at 
8:00 UT on November 20. Interestingly, the northern pore faded 
away in parallel to this growths spurt, (i) Most of the penumbra 
in the leading spot has disappeared by 5:00 UT on November 21. 
(j) Finally, on November 23, all spots, pores, and magnetic knots 
had vanished and only a bright plage region remained until it ro- 
tated off the visible hemisphere. 

In summary and neglecting all details of active region evolu- 
tion, active region decay rates can be computed using a linear fit 



for time periods when the area coverage reached its maximum to 
the point when the area fell below 5 Mm^. The overall decay rate 
of the active region is 72.6 Mm^ per day. The values for penum- 
brae and umbral cores/pores are 48.3 and 25.1 Mm^ per day, 
respectively. Similarly, we computed the growth rate of umbrae 
and penumbrae starting at 21:00 UT on November 14. The val- 
ues are 171.4 and 104.8 Mm^ per day, respectively. This is about 
four times faster than the corresponding decay times but less 
than one half of the penumbral growth rate of about 400 Mm^ 
per day presented by Schlichenmaier et al. (2010). 

Hinode G-band images allow us to zoom in on the two trail- 
ing spots. The region-of-interest (ROI) is shown in the top row 
of Fig. 5. The data of 2010 November 18-22 are corrected for 
geometrical foreshortening and the center-to-limb variation. The 
ROI with a size of 756 x 756 pixels or 60 Mm x 60 Mm was 
centered at a heliographic latitude of S3 1.9°. The first G-band 
image of the daily observing sequences was used as a reference 
to align all other data. 

Various solar features are identified using intensity thresh- 
olds and morphological image processing. We indiscriminately 
used a fixed intensity threshold of /n,ag = 0.8/() for strong mag- 
netic features and an adaptive threshold for G-band bright points 
of /bp - (1.37 - 0.08/i)/o, where fi is the cosine of the helio- 
centric angle 9, and Iq refers to the average quiet Sun intensity. 
Intensity values between /bp and /^ag consequently enclose gran- 
ulation. The measured spot areas for HMI and Hinode agree with 
each other, and the remaining differences can be attributed to dif- 
ferent spectral characteristics of the observed passbands, image 
scales, and threshold/selection criteria. In general, the temporal 
evolution within the ROI follows the same trend as discussed 
in the context of SDO observations (Fig. 3). Areas, velocities, 
and other physical quantities based on high-resolution data are 
included in Tab. 1 . 

In the following, we present a chronology of the important 
phases of sunspot decay based on the high-spatial resolution 
G-band images. On November 18, the two trailing sunspots 
were embedded in a network of G-band bright points hinting at 
widely dispersed, weak magnetic fields. Both spots were filled 
with numerous umbral dots. Spot A had three umbral cores sep- 
arated by faint light bridges and a small penumbral segment 
pointing westward. Spot B had two umbral cores with a few 
associated magnetic knots. It possessed curved penumbral fil- 
aments pointing eastwards, which are indicative of twisted mag- 
netic field lines. These kind of non-radial penumbral filaments 
are frequently observed in flaring sunspots with horizontal shear 
flows (Denker et al., 2007; Deng et al., 2006). On November 19, 
the faint light-bridges had disappeared and the penumbra had 
vanished leaving only a single pore, which was filled with con- 
glomerates of umbral dots. Spot B also lost its penumbra leav- 
ing four umbral cores separated by faint light-bridges, which, 
however, were more pronounced in comparison to previous day. 
Magnetic knots were still surrounding spot B. Spot A had de- 
cayed on November 20 with only two faint magnetic knots re- 
maining at its location. Interestingly, at this point of time spot B 
started to grow in area with hints of penumbral filament being 
visible on its western side. It consisted of three umbral cores, 
which were now separated by strong light-bridges, which split 
the spot in two halves along its north-south axis. The presence 
of strong light-bridges might suggest the initiation of the spot's 
decay phase (see Sobotka et al., 1993). Once spot A had disap- 
peared on November 21, it did not leave any significant trace 
within the network of G-band bright points. By this time spot B 
had dissolved into multiple tiny pores and magnetic knots, which 
roughly covered the same region as on the previous day. By 
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Fig. 5. G-band images {top) after correction of geometrical foreshortening tracing the photospheric evolution of the central part 
of active region NOAA 11126 {from left to right 2010 November 18-22). The horizontal flow speeds {2nd row) are given by the 
velocity scale to the right. The direction of the horizontal flows {3rd row) are displayed according to the color-wheel and arrows 
for which a velocity of 1 km s ' corresponds to exactly the grid spacing. The divergence of the horizontal flow field {bottom) are 
presented according to the scale on the right, where gray indicates divergence values close to zero. All G-band images and LCT 
maps were aligned so that the center of the panels coincides with a latitude of 3 1 .9° South. The white and black contours outline the 
location of the small sunspot/pores. Times refer to the first image of the time-series, which was used to compute the flow maps. 



November 22, only two tiny pores were left from spot B, which 
disappeared on November 23 just leaving G-band bright points 
at its point of disappearance. Two observations are noteworthy: 

(1) The area covered by the G-band bright points remained al- 
most constant during the disk passage of the active region, which 
suggests that flux decays more or less in place and is not redis- 
tributed over a larger area. The time scale of flux removal or dis- 
persal extends well beyond the photometric decay time of strong 
magnetic features such as sunspots, pores, and magnetic knots. 

(2) The two trailing spots have different histories of flux emer- 
gence and decay. 

3.2. Chromospheric evolution 

The description of the chromospheric evolution is based on 
CaiiH (top row of Fig. 6) and Ha line core intensity maps (top 
row of Fig. 8). The CaiiH images were averaged over three 
hours to highlight some of the long-lived chromospheric fea- 
tures. On November 18, the lower quarter of the CaiiH image 



showed the undisturbed pattern of inverse granulation (Rutten 
et al., 2004). Ca ii H brightenings cover a larger area, since their 
filling factor is significantly greater than that of G-band bright 
points. Spot A was encircled by individual brightenings at a ra- 
dial distance of 5 Mm. These individual brightenings coalesced 
into a wagon-wheel-like Ca ii H intensity structure in the three- 
hour average image pointing to the presence of MMFs. A similar 
feature was noticeable around a much smaller pore with a diame- 
ter of about 2 Mm located to the west of spot A. This conspicuous 
Ca II intensity structure was first described by Shine et al. ( 1 996), 
who identified the bright ridges with azimuthal convergence re- 
gions, whereas the dark regions between the spokes correspond 
to azimuthal divergence zones. Spot A resides in the middle of 
a supergranular cell with a diameter of about 20 Mm. In con- 
trast, spot B is embedded in an area of much more pronounced 
Ca II H brightenings. The exterior tips of its penumbral filaments 
appeared bright because a low-level B-class flare occurred at the 
time of observations. 
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Fig. 6. CaiiH average images tracing the evolution of active region NOAA 11126 in the chromosphere. The horizontal proper 
motions were derived from time-sequences of CaiiH images. Otherwise, data processing and display are the same as in Fig. 5. 
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Fig. 7. Radial dependence of parameters char- 
acterizing the flow fields of spots A (left) and 
B (right). Radial distances are measured from 
the outer boundary of the spots, i.e., 0.0 Mm 
marks the transition from penumbra to gran- 
ulation. The thick black curves indicate the 
average angle of horizontal flow vectors with 
respect to the outward radial direction. High 
values refer to inward flows and low values 
indicate outward flow vectors. Fractions of in- 
/outward flow vectors as defined in the text 
are shown as a dash-dotted and dashed lines, 
respectively. The horizontal straight lines de- 
notes the 100%-level. The average flow speed 
is displayed as a dash-dot-dot-dotted line. 



On November 19, the average CaiiH image showed evi- 
dence that the moat flow around spot A survived the initial 
stages of sunspot decay, even after its penumbra had disap- 
peared (see e.g., Deng et al., 2007). Spot A had completely dis- 
solved by November 20. Assuming that spot A was located in 
the center of a supergranular cell as indicated by the surround- 
ing CaiiH brightenings on November 18 and 19, we conclude 
that this supergranular cell ceased to exist on November 20, 
when the strong magnetic fields of spot A where no longer 
present. The Ca n H brightenings are now squeezed together by 
two supergranular cells to east and to the west, which could al- 
ready be identified on November 19 and remained visible until 
November 22. 

At the time when spot A vanished, spot B had already frag- 
mented into two umbral cores, which were separated by a strong 
light-bridge with a noticeable dark lane along its axis (Rouppe 
van der Voort et al., 2010). This dark core was even more clearly 
discernible in the three-hour average CaiiH image. The pres- 
ence of strong light-bridge can be taken as a first indication of 
the spot entering the decay phase. During the next two days on 
November 21-22, spot B decayed further 

The Ho- line core intensity map of November 18 exposed 
a radial pattern of fibrils around spot A reminiscent of a super- 
penumbra. Strong brightenings in the Hor line core in the vicinity 



of spot B are related to the aforementioned small B-class flare. 
We did not notice any significant Yia intensity features on the 
following days, except that after the decay of both spots the Ha 
plage became more prominent. By the end of our observations 
with the Echelle spectrograph, only a very compact plage region 
with a length of about 8 Mm remained within ROI. 

3.3. Horizontal proper motions 

The LCT flow maps displayed in Figs. 5 and 6 were computed 
using the time-sequence of Hinode G-band and CaiiH images. 
As mentioned in Section 3.1, we distinguished various solar fea- 
tures (e.g., G-band bright points, granulation, strong magnetic 
features, and sunspot penumbrae) by morphological and adap- 
tive thresholding. For comparison, we refer to the typical flow 
speeds of granulation Vgran = 0.47 + 0.27 km s"' and G-band 
bright points Vbp = 0.23 + 0.15 km s"' as reported by Verma & 
Denker (201 1), who also presented values for longer time inter- 
vals Ar (Tab. 1 therein) and cadences Af (Fig. 3 therein). Their 
values for granulation are in very good agreement with the mean 
values Vgran = 0.36 + 0.20 km s"' of the present study. The mean 
values for bright points Vbp - 0.26 + 0.15 km s"' in the neigh- 
borhood of active region NOAA 1 1 126 are within the range of 
expected values. Daily values of flow speeds for granulation and 
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Fig. 8. VTT Echelle spectra were used to study the chromospheric evolution of active region NOAA 111 26 during the period from 
2010 November 18-22. The chromospheric fibril structure surrounding the decaying sunspots can be traced in Ha /i656.28 nm 
line core intensity {top) and LOS velocity {middle) maps. In addition, photospheric LOS velocities {bottom) were derived from 
Hinode/SP Fe i /1630.25 nm spectra. The FOV matches those of Figs. 5 and 6. Regions not covered by the spectral data are shown 
in medium gray. The times above the panels refers to the start of a spectral scan. 



G-band bright points are given in Tab. 1. The standard devia- 
tion of the aforementioned horizontal flow speeds refers to the 
variance in the data for a particular solar feature rather than to a 
formal error The intrinsic error of the LCT algorithm is 50 m s ' 
and 15° for flow speed and direction, respectively (see Verma & 
Denker, 2011). 

We computed the flow vectors for strong magnetic elements, 
i.e., spots A and B. We identified spots by smoothing the ge- 
ometrically corrected G-band images using a Gaussian kernel 
with a FWHM of 1280 km and subsequently applied an inten- 
sity threshold of 0.83/o. We plotted contours based on these 
binary masks in Figs. 5 and 6 to provide some visual guid- 
ance in identifying the spot positions in the physical maps. The 
mean flow speeds vspotA - 0.19 ± 0.17 km s"' and vspot b = 
0.18 + 0.11 km s"' are virtually identical for both spots. While 
computing the mean velocity v>spotB for spot B, we discarded 
November 22, because on that day only a tiny pore with a di- 
ameter of 4 Mm was left at the location of spot B. Therefore, ad- 
jacent regions with G-band bright points and granulation would 
bias the flow speed towards higher values. The overall flow pat- 
tern for the Ca ii H data is very similar as compared to G-band 
data, and the average values vspot a - 0.18 + 0.20 km s"' and 
vspotB - 0.18 + 0.08 km s"' are virtually identical with a ten- 
dency to be slightly lower on individual days. 

The moat flow around sunspots and pores reveals itself as 
radially, outward-directed vectors in flow maps, which point to 



a ring-like structure with kernels of elevated flow speeds. In the 
speed and azimuth maps on November 18 and 19 (Figs. 5 and 6) 
indications of moat flow were visible around spot A and a small 
pore located to the west of spot A. There was no clear signature 
of the moat flow detectable around spot B. To express the charac- 
teristics of the moat flow in more quantitative terms, we plotted 
in Fig. 7 the angle of the flow vectors, the fraction of in- and 
outward flow vectors, and the horizontal velocities for spots A 
and B. These quantities are radial averages and the zero point 
of the radial distance was placed at the penumbra/granulation 
boundary of the spots. We considered a flow vector to point in- 
/outward, if the tip of the arrow aims in-/outward and if the angle 
with the radial direction is less than 20°. Therefore, the frac- 
tions of in-/outward flow vectors presented in Fig. 7 do not add 
up to 100%, since more tangential flows are neglected. At a ra- 
dial distance of 6 ± 2 Mm from the boundary of spot A more 
than 50% of flow vectors point outward. The average flow speed 
in this region is about 0.5 km s"'. The region with high-speed 
outward flows marks the location, where the moat flow around 
spot A terminates. An exact determination of where the moat 
flow starts and where it ends strongly depends on the underly- 
ing criteria, e.g., photometric or magnetic features, or horizontal 
proper motions. In Sect. 3.5, we will discuss some properties of 
moat flow based on the proper motions of MMFs. Sobotka & 
Roudier (2007) determine that the moat radius is independent of 
the spot size. However, Balthasar & Muglach (2010) find that the 
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moat flow terminates at a distance of four times the spot radius 
- in contrast to three times the spot radius in the present study. 
Extended statistical studies will help to clarify this issue. 

We created two ring-like structures by morphological dila- 
tion with a width of 4 Mm, which encircled both spots at a dis- 
tance of 2 Mm. This ring-like structure starts, where inflows turn 
to outflows 2 Mm beyond the boundary of spot A, and termi- 
nates, where the outflows reach the highest speed. We used these 
templates to calculate the flows in the immediate neighborhood 
of spots A and B. These regions were labeled as rings A and B. 
The velocities for these regions are included in Tab. 1 . The flow 
speed VRing a was more than double than that of vspot a because 
of the suppressed convective motions in the spot's interior. The 
flow speeds in the immediate vicinity of spot B was higher by 
30-40% on November 18 and 19. However, on November 22 
during the final decay stage of spot B, no major difference be- 
tween the spot and its closest surroundings were observed. 

To effectively visualize the plasma motions we presented 
high resolution maps of flow speed and azimuth (see Figs. 5 
and 6). The speed maps of G-band and CaiiH images are vir- 
tually same. The region with G-band bright points exhibits sup- 
pressed velocities in both cases. The regions around spot A and 
the small pore to the west are surrounded by a ring of high ve- 
locities (around 0.5 km s '), where the ordered moat flow ter- 
minates. In the azimuth maps of November 18 and 19 outward 
plasma motions are traceable in these regions. In addition, to 
gain more insight into inflows and outflows around the spots we 
computed the divergence of the flow field. The divergence maps 
are included in the bottom row of Fig. 5. The positive values of 
divergence refer to outflows, whereas negative values indicate 
inflows. Negative divergence values are encountered within the 
boundaries of the sunspots. In contrast small patches of positive 
divergence encircle the spots. This is indicative of inflows in the 
sunspots and outflows at their periphery. The area with low di- 
vergence values is increasing as the active region is decaying. 
The presence of a low divergence region could be a signature of 
the final stages of decay. 

3.4. Line-of-sight velocities 

The LOS velocities were derived using Hinode/SP Fe i spectral 
data and Ha spectra of the VTT Echelle spectrograph. LOS ve- 
locities for the Fe i spectral line were calculated using the Fourier 
phase method (Schmidt et al., 1999). This method makes use 
of entire line profile and is less sensitive to noise, along with 
this it takes into account the spectral line asymmetry. To com- 
pute the LOS velocity for Ha spectra we calculated shifts us- 
ing parabola fits to the central 50 pixels (0.03 nm), because the 
Ha spectral line is too wide to identify the real continuum. The 
calculated shifts in both cases were converted to velocities us- 
ing the Doppler formula. The average photospheric velocity of 
dark umbral cores was used as the frame of reference. The Fe i 
/1656.92 nm line served as the reference for the Ha velocities. 
The Doppler velocity maps are displayed in the bottom and mid- 
dle rows of Fig. 8. Redshifts in these maps are positive and 
blueshifts are negative, hence, areas moving away from the ob- 
server are bright, while areas moving toward the observer are 
dark. 

To compute velocities for various solar features we used an 
intensity mask generated using G-band images. The values are 
compiled in Tab. 1. On November 18 and 20 in the Fei Doppler 
velocity maps strong photospheric downflows were observed at 
edges of spot B. The average downflow velocity in these regions 
is about 2.5 km s"' . In all maps a gray patch of near zero veloc- 



ity was observed in the central FOV, which became more promi- 
nent on the last two days of observation. We identify this region 
with abnormal granulation (see e.g., de Boer & Kneer, 1992), 
in which convection is still strongly inhibited by the presence of 
(dispersed) magnetic fields. 

In case of Ha LOS velocities no conspicuous features were 
visible in the velocity maps, except on November 18 when we 
observed a radial filamentary structure around spot A in Ha 
line core intensity maps, which resembled a superpenumbra. At 
the footpoints of the dark Ha filaments, we measure downflows 
from about 3.5 up to 4.5 km s"' , which we interpret as an inverse 
Evershed flow (Maltby, 1975). Tightly wound superpenumbral 
spirals are only predicted for spots with radii larger than 8 Mm, 
whereas spot A was compact with a radius of 4 Mm. The down- 
ward chromospheric velocities at the edge of spot A are compat- 
ible with the MHD model of superpenumbral flows presented by 
Peter (1996). 

3.5. Magnetic fields 

In addition to the photometric evolution shown in Fig. 3, we 
also computed the flux contained in the active region as shown 
in Fig. 4. Since only HMI LOS magnetograms were available 
(processing of the vector magnetograms is still under way). We 
took the measured magnetic field strength at face value and 
only carried out a geometrical correction to yield the proper av- 
erage values of the magnetic flux. The geometrical correction 
only applies to the surface area covered by a pixel, which can 
simply be achieved by dividing the magnetic field strenght by 
fi = cos6. Signatures of geometric projection eff'ects can be seen 
when the active region was close to the east limb, and the an- 
gle between LOS and shallow penumbral field lines leads to an 
apparent polarity reversal. Consequently, the positive and neg- 
ative flux gradients have opposite signs, while to total flux re- 
mained almost constant for the first 30 hours until projection 
effects become less severe. To compute the temporal evolution 
of magnetic flux, we created a binary template, which only con- 
tained pixels above/below +50 Gauss in the HMI magnetograms. 
Morphological erosion with a 1-Mm kernel was applied to the 
template three times to eliminate small, isolated flux concen- 
trations. Finally, we used morphological dilation with a 5-Mm 
kernel to include the strong magnetic fields in the immediate 
neighborhood of active region NOAA 11126. In this way, we 
avoided a bias, which could be introduced by weak magnetic 
fields, which do not belong to the active region but are contained 
within the FOV shown in Fig 2. 

The growth rates of the magnetic flux are 2.66, 1.98, and 
1.09 X lO'-' Wb day"' for the total, positive, and negative flux, 
respectively. The negative flux showed a monotonous rise un- 
til about 06:00 UT on November 16, whereas the positive flux 
increased with a twice steeper slope, stopped a day earlier, and 
turned to a shallower slope until the end of November 17. On 
average the positive flux was three times larger than the negative 
one. However, using adaptive thresholding and morphological 
image processing tools, we only measured the flux in the imme- 
diate vicinity of the sunspots. The missing negative flux required 
for flux balance has to be contained in flux concentration beyond 
the immediate neighborhood of the sunspots. The decay rates of 
the magnetic flux are four to five times lower than the growth 
rates and amount to 0.66, 0.47, and 0.23 x 10'^ Wb day"' for the 
total, positive, and negative flux, respectively. A linear fit to the 
data was used and despite some deviations from a linear trend, 
there is no indication for a parabolic (or any other) decay law. 
The decay rates of the present study are in agreement with Kubo 
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Fig. 9. Maps of physical parameters derived using SIR-code for Hinode/SP scan on 2010 November 18 (from top-left to bottom- 
right) normalized intensity ///q, vertical component of magnetic flux density B^., magnetic field azimuth (f), Doppler velocity vlos, 
horizontal component of magnetic flux density Z?hoi , and magnetic field inclination j. 



et al. (2008) who report rates on the order of 10'^ Wb day ' and 
discusses the magnetic flux loss rate in a decaying active region. 

The SIR-code (Stokes Inversion based on Response func- 
tions) developed by Ruiz Cobo & del Toro Iniesta (1992) was 
used to invert the Hinode/SP spectra. We restrict ourselves to 
the more sensitive line Fei /I630.25nm (geff - 2.5). The start- 
ing model covers the optical depth range -i-l.O < logr < -4.4. 
A limb-darkening factor is considered according to Eq. 10 of 
Pierce & Slaughter (1977). We assume a constant macroturbu- 
lence of 1 kms"' and a fixed stray light contribution of one per- 
cent. The inversions deliver the temperature stratification with 
four nodes T{t), the total magnetic flux density Btot, the mag- 
netic inclination y and azimuth 0, and the Doppler velocity vlos 
constant with height (one node for each of these physical param- 
eters). 

The magnetic azimuth ambiguity must be solved after the 
inversions. For the first two days when we observed two ma- 
jor pores/small sunspots, we assumed two azimuth centers (see 
Balthasar, 2006). For the other days it is sufficient to assume 
one azimuth center away from the pores in the direction towards 
disk center. The magnetic field in such small pores is more or 
less vertical to the solar surface so that on these days the field 
is sufficiently inclined with respect to the LOS that such an as- 
sumption is justified. If the expected azimuth deviates by more 
than 90° and less than 270° from the calculated one, we correct 
it by 180°. Finally, we rotate the magnetic vector with respect to 
the local solar frame. For a few locations, it happens that there 
is a sudden change of sign in the Cartesian components of the 
magnetic vector This problem is solved by an additional correc- 
tion of the LOS-azimuth. To be on the safe side, we additionally 



apply to data of November 1 8 the code provided by Leka et al. 
(2009), which minimizes + |VB| to solve the azimuth ambi- 
guity. Pixels, where the integrated circular and the integrated lin- 
ear polarization are below 0.006, are excluded from the further 
analysis of the magnetic field. We finally use the magnetic vec- 
tor field in the local solar frame. The results of the SlR-inversion 
for November 18 are depicted in Fig. 9. 

The vertical magnetic flux density B~ in the ROI is predomi- 
nantly positive and points outwards. The only significant patches 
of negative polarity were found on November 18 on the east- 
ern side of spot B. Penumbral filaments connect spot B to sev- 
eral magnetic knots of negative polarity. In addition, we measure 
Evershed flows with velocities vlos > 2 km s"' . On smaller spa- 
tial scales, we find MMFs in the vicinity of spot A. They are 
mostly of type 11 (unipolar with the same polarity as the spot) 
but a few scattered U-shaped type 1 MMFs (bipolar with the in- 
ner footpoint of opposite polarity of the spot) were observed as 
well. No type 1 MMFs were observed near spot B. 

Since the time cadence of Hinode/SP scans is about 12 min, 
we used SDO/HMl magnetograms with a cadence of 45 s to 
study the dynamics of MMFs. Space-time slices are an option 
to visualize local changes of the magnetic flux in a time-series. 
The continuous coverage of HMl magnetograms allowed us to 
depict in Fig. 10 the temporal evolution of the magnetic flux 
around spot A during a 12-hour period on 2010 November 18. 
The magnetic flux changes were recorded within two 0.5-Mm 
wide annuli of 5.75 Mm (bottom panel) and 7.25 nm (top panel) 
radius, respectively. Note that the linear scale of the the ordinate 
in Fig. 10 is 38 Mm and 47 Mm for the inner and outer ring, 
respectively. The inner ring corresponds to the location, where 
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Fig. 10. Space-time slices showing the tem- 
poral evolution of the magnetic flux along 
annuli with radii of 5.75 Mm {bottom) and 
7.25 Mm (top), respectively, which are cen- 
tered on spot A. HMI magnetograms with a 
cadence of 45 s covered a period of 12 hours 
starting at 04:00 UT on 2010 November 18. 
The annulus with a widths of 0.5 Mm was 
cut open in the East and then transformed to 
straight line so that the southern part of the 
spot is mapped to the lower half of the space- 
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Fig. 11. Local, background-subtracted variation of the magnetic 
flux Bvar around spot A on 2010 November 18. 



the strongest outflows were observed (see Fig. 7). The outer ring 
marks the site, where the largest fraction (about 8%) of the an- 
nulus is covered by negative-polarity features. 

Within a radial distance from 4-10 Mm from the center of 
spot A, only about 5% of the magnetic flux observed during 12- 
hour period is of negative polarity. Negative-polarity MMFs are 
weaker by about 20% as compared to MMFs of opposite polar- 
ity. Negative flux elements are most prominent in the southern 
half of spot A and absent in the northern half. The most notable 
difference between positive- and negative-polarity features is the 
morphology. Negative-polarity features can be tracked in space- 
time slices for only about 15-30 min, i.e., they are strongly local- 
ized. Time-lapse movies show that these negative-polarity fea- 
tures are mostly type I MMFs. The positive polarity appears first 
in the space-time slices (to the left) followed by the negative 
polarity. The appearance of the positive-polarity features is very 
different. Horizontal striation can be seen for several hours. They 
can abruptly (dis-)appear and in a few cases merging as well 
as splitting of the striae can be observed. In time-lapse movies, 
type II MMFs are moving away from spot A along radial paths 
on top of a background of positive flux. 

The temporal variation of the magnetic flux above/below the 
local background is given by Byar = (1^ - {B)\), where (. . .) in- 
dicates a temporal average. To compute the temporal average in 



Fig. 11, we used the same 12-hour time-series of HMI magne- 
tograms that was used to create Fig. 10. The strongest changes 
occur at the periphery of spot A, where also the horizontal mag- 
netic field Bhor is strongest as can be seen in Fig. 9. This corre- 
spondence also applies to the feature of strong horizontal field to 
the north-east of spot A, which belongs to the (chromospheric) 
network. The most notable features in Fig. 1 1 are, however, sev- 
eral spoke-like structures with a length of about 10 Mm, which 
extend from spot A all the way out to the border of the supergran- 
ule surrounding the spot. These elongated structures are the sig- 
natures of MMFs streaming from the spot along identical paths 
(Harvey & Harvey, 1973) towards the supergranular bound- 
ary. However, the distance traveled is much shorter as Harvey 
& Harvey (1973) observed for mature sunspots. Therefore, we 
conclude that the preferred paths of MMFs are a phenomenon, 
which can be observed even in the final stages of sunspot de- 
cay. The traces left by MMFs in the variance of the magnetic 
field around its background Bvar are most prominent on the east- 
ern side of the spot. They also match the wagon-wheel struc- 
ture discussed in the context of the time-averaged Ca ii H images 
in Sect. 3.2. Finally, the preference of MMFs for distinct radial 
channels connecting the sunspots' magnetic fields to the strong 
network fields was also noticed by Hagenaar & Shine (2005). 

Maps similar to Fig. 9 were computed for each day so that 
we could study the evolution of the magnetic field. The strongest 
horizontal magnetic flux densities Z?hoi are observed in the im- 
mediate surroundings of spots A and B as long as the spots are 
compact (November 18 and 19). They symmetrically enclose 
the entire boundary of the spots. Once spot A had disappeared 
and spot B started to fragment, this symmetry is broken most 
evidently on November 20, when strong horizontal fields were 
present at the western side of spot B. This corresponds to a time 
when a rudimentary penumbra was present. Whenever indica- 
tions of penumbral filaments were observed during the decay of 
the sunspots, inclined magnetic field lines and Evershed flows 
were present at the same time (cf., Leka & Skumanich, 1998; 
Yang et al., 2003). 

The magnetic field lines spread out symmetrically from the 
center of the spots on November 18 and 19. Such well defined 
azimuth centers can still be found in the later decay stages of 
spot B. Even though these region are no longer circular and be- 
come elongated. Smaller azimuth centers are also observed to 
the east of spot B on November 19. These centers are associated 
to several magnetic knots. 

The Doppler velocity vlos is suppressed in the presence of 
strong magnetic features. Starting on November 20, the veloc- 
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ity pattern associated with the magnetic fields changes. Positive 
Doppler velocities occupy more and more of the magnetic re- 
gion. This coherent pattern was also observed in the divergence 
maps (bottom row of Fig. 5) and the horizontal flow maps 
(Figs. 5 and 6). Indeed, considering that the region is approach- 
ing the limb, some of the Doppler velocities can be interpreted 
as a coherent proper motion of the magnetic region towards the 
south-west (cf. the azimuth maps of the horizontal proper mo- 
tions). In summary, during the final stages of sunspot decay, the 
three-dimensional flow field, in regions previously occupied by 
strong magnetic fields, significantly differs from granular flow 
patterns or regions of reduced velocities in the presence of G- 
band bright points. Since we are only presenting a case study, 
the question remains, if this intriguing flow pattern is a typical 
feature of sunspot decay. 

4. Conclusions 

We have presented a detailed account of the final stages in the 
decay of the active region NOAA 1 1 126, which did not obey the 
Hale-Nicholson polarity law (e.g., Zirin, 1988). Since only one 
out of ten active regions shows such a behavior (Howard, 1990) 
and we only present a case study, our results might not be repre- 
sentative for sunspot decay in general. However, space missions 
such as Hinode and SDO provide nowadays data of sufficient 
coverage, resolution, and cadence that statistical properties of 
sunspot decay become accessible. Furthermore, previous studies 
of non-Hale regions (e.g., Lopez Fuentes et al., 2000, and ref- 
erence therein) were centered on flux emergence, (5-spots, and 
strong solar flares. The present study can consequently be con- 
sidered as an extension of these studies with a focus on a much 
quieter magnetic field topology, which might be representative 
for the lower solar activity during cycle No. 24 (Petrovay, 2010; 
Nielsen & Kjeldsen, 2011). 

The major findings of our study can be summarized as fol- 
lows: (1) MMFs were observed in the vicinity of spot A un- 
til it decayed. Mostly type II and a few interspersed U-shaped 
type I MMFs contributed to the observed moat flow, which 
also left a clear signature in the time-averaged CaiiH images 
(Martinez Fillet, 2002). (2) Even though penumbral filaments 
had almost completely disappeared in photospheric G-band im- 
ages of spot A on November 18, \ia line core images clearly 
exhibited a structure reminiscent of a superpenumbra. Thus, fila- 
mentary structures including the inverse Evershed flow (Maltby, 
1975; Georgakilas & Christopoulou, 2003) might be more per- 
sistent in the chromosphere. (3) We have also observed MMFs in 
the vicinity of a tiny pore with a diameter of about 2 Mm, which 
did not show any indication of penumbral filaments. Such an ob- 
servation argues strongly against a close tie between Evershed 
flows and MMFs (cf., Vargas Dommguez et al., 2008, 2010; 
Cabrera Solana et al., 2006). Rempel (2011) argues based on 
MHD simulations of sunspots that penumbral flows can be sepa- 
rated in two components, where the shallow one corresponds to 
the Evershed flow and the deeper one is related to moat flow. (4) 
The strong rotation and twist seen in spot B might explain, why 
this trailing spot never advected sufficient magnetic flux to es- 
tablish more than a rudimentary penumbra and remained highly 
fragmented during its entire life cycle. (5) The photospheric and 
chromospheric maps of horizontal flows show a peculiar pattern, 
once the last dark feature of the active region has disappeared. 
In general, the flow field in this region is less structured than re- 
gions covered by granulation, i.e., the dispersed magnetic field 
still significantly affects the convective pattern. Chromospheric 
flows have increased notably compared to times when spots and 



pores were still present. Most prominently, a contiguous area of 
low divergence appears towards the end of sunspot decay. 
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